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1. Introduction 

Long-range interactions are commonly defined as pair interactions decaying fike r~" , where r is the distance 
between the interacting constituents, and the positive exponent a is smaller than or equal to the spatial dimension 
d of the system. Examples of long-range interactions are common in nature, including gravitating masses, or 
Coulomb forces in the absence of screening. It has been known for several decades that long-range systems may 
show unconventional and sometimes surprising behaviour compared to the more familiar short-range interacting 
ones. For example, negative heat capacities were shown to exist in bounded self-gravitating gas spheres by 
Lynden-Bell and Wood [1] in 1968, and Thirring explained this observation by relating it to the nonequivalence 
of microcanonical and canonical ensembles in the presence of long-range interactions [2]. It was only recently in 
[3, 4] that nonequivalence of ensemble has also been observed (theoretically) in long-range quantum spin systems 
like the one studied in the present article. Another particularly interesting phenomenon, going under the name of 
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quasistationary states, has triggered intense activity, reviewed in [5], in the field of long-range interacting systems. 
The term quasistationary is used to describe metastable states whose lifetime diverges with increasing system 
size A'^. The physical importance of quasistationary behaviour should be obvious: For a sufficiently large system, 
the transition from a quasistationary state to equilibrium takes place on a time scale that is larger than the 
experimentally accessible observation time. Hence, equilibrium properties will not be observed, and instead the 
statistical properties of the quasistationary regime are of interest. 

Most of the studies of quasistationary states have focused on the Hamiltonian Mean-Field model [6], a toy model 
consisting of classical XY spins (or plane rotators), each coupled to every other at equal strength. The typical 
scenario observed is that, for a suitable class of initial conditions, the total magnetization of this spin model 
rapidly relaxes to some quasistationary value different from its equilibrium value, remains there for a very long 
time, and finally decays to equilibrium (sec Figure 33 in [5]). More recently, quasistationary states have been 
observed in gravitational sheet models [7] and their existence has been argued to be generic for a large class of 
classical long-range systems [8]. Virtually all finite-system results in the field have been obtained by numerical 
techniques, supplemented by analytical calculations in the TV — > oo Vlasov continuum limit [5, 9-12]. 
In a recent Letter [13] , the author has reported results on the dynamics of a long-range interacting version of the 
Emch-Radin model [14, 15]. This quantum mechanical spin model is particularly amenable to analytic calculations 
and is known to show relaxation to equilibrium in a non-Markovian way [14] . Unlike in the papers of Emch and 
Radin, long-range interaction are permitted, i.e. exponents a (determining the spatial decay of the interaction 
strength) which are smaller than the lattice dimension d. The results announced in [13] include an analytic 
expression for the expectation value {A){t) as a function of time t for a certain class of observables A, as well as 
upper and lower bounds on {A){t) in the thermodynamic limit of infinite system size. In particular, the lower 
bound establishes rigorously the occurrence of quasi-stationary behaviour in the long-range version of this model. 
The aim of the present article is twofold: (a) to provide a detailed derivation of the bounds on {A){t), and (b) to 
numerically investigate how the equilibration time scales with the system size A'^. 

The article is structured as follows: In Section 2, the Emch-Radin model is introduced, generalizing also to long- 
range spin-spin interactions. In Section 3, an analytic expression for the time evolution of the expectation value 
{A){t) is given for a class of observables A and initial states p(0). This is a rather straightforward extension of 
Emch's analysis to the long-range case. Section 4 is devoted to the study of {A)(t) in the thermodynamic limit of 
infinite system size. In Section 4.1, an upper bound on {A){t) is constructed. For exponents < a < 1, the lower 
bound derived in Section 4.2 establishes that, in the thermodynamic limit, the system is trapped for all times at 
its initial state and relaxation to equilibrium is never observed. In Section 4.3 we perform a numerical analysis of 
the exact formula for {A){t), studying in particular the scaling behaviour of the equilibration time scale To with 
the system size. We find a power law To oc A*"* with some exponent q, in agreement with what is predicted by the 
upper bound derived in Section 4.1. The findings are summarized and discussed in Section 5. 
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2. Emch-Radin model 

Consider N identical spin-1/2 particles, attached to the sites {1, . . . , A'"} of a finite one-dimensional lattice with 
periodic boundary conditions.^ The quantum dynamics takes place on the Hilbert space 

JV 

n = (g)cl (1) 

i = l 

where the Cf are identical replicas of the two-dimensional Hilbert space of a single spin-1/2 particle. The time 
evolution on H is generated by the Hamiltonian 

JV N/2 N 

Hn =JVa'^Yl -h^ol, (2) 

i = l j = l 1 = 1 

with (j^ denoting the z component of the Pauli spin operator and /i £ R an external magnetic field in the z 
direction. The index i -\- j is to be considered modulo TV to account for periodic boundary conditions. The e[j) 
are pair coupling constants depending on the distance j of two spin operators on the lattice, and we assume 

lim e{j) = 0. (3) 

A typical example we have in mind is the Dyson model [16] with algebraically decaying couplings, e(j) = j~" , 
but, in contrast to Dyson's work, exponents a smaller than 1 are also allowed. The explicit calculations reported 
in the present article are for such Dyson- type interactions, but generalizations are possible. For exponents a > 1, 
the interaction is absolutely summable, 

I]k(i)l<«=, (4) 

and this case has already been treated in [14, 15, 17]. Here we are interested in long-range interactions with 
< a < 1 where the sum in (4) diverges. In this case, the double sum in (2) gives rise to an infinite energy per 
spin {Hn)/N in the thermodynamic limit A'^ — oo. To render this limit finite, the normalization 

/ N/2 \ -1 

^^=[2j2<j)j (5) 

has been introduced in (2). This normalization is a generalization of the so-called Kac prescription (introduced 
by Baker [18]) commonly used when studying long-range interacting systems. For large system sizes A^, its 



^ Generalizations to higher- dimensional lattices are straightforward, but come at the expense of a more complicated 
notation. All results reported in this article generalize to dimensions d > 2 with only insignificant modifications. 
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asymptotic behaviour is of the form 



■q: — 1 



for < a < 1, 



2A/'c« ~ < l/lnAT 



for a = 1, 



(6) 



for a > 1, 



where ^ denotes the Riemann zeta function. Equation (6) implies that J^a vanishes in the thermodynamic limit 

for exponents a < 1. 

We consider observables of the type 



with real coefficients a^, i.e. observables which are linear combinations of the x components erf of the spin operators. 
None of these observables commutes with the Hamiltonian (2), and hence there is a chance of observing relaxation 
to equilibrium when studying the time evolution of expectation values {A}{t). For technical reasons, the analysis 
is restricted to initial states (i.e. initial density operators) p(0) which are diagonal in the erf tensor product 
eigenbasis of H. This assumption simplifies the calculations but is not expected to be essential for the findings 
reported. 

The Emch-Radin model, including the choices of observables and initial states, is inspired by induction decay 
experiments probing the pulsed magnetic resonance of nuclei in a CaF2 crystal [19]. In these experiments, the 
decay of the x component of the total magnetization ^(1, . . . , 1) = X^ili ""f measured, and it was found 
to be superimposed by oscillations (see Fig. 1 of [19]). Despite its simplifying assumptions, the model captures 
well both the decay and the beating observed experimentally. Moreover, the Emch-Radin model and its random 
generalizations [17] have proved useful as paradigmatic models for which the approach to equilibrium in quantum 
systems can be studied analytically. 

3. Finite-system time evolution 

We study the time evolution, generated by the Hamiltonian Hn, of the expectation value of an observable A of 
the form (7) with respect to the initial state p(0). 



(in units where H = 1). Performing the trace in the erf eigenbasis of "H, the diagonal form of p(0) implies that only 
the diagonal elements of the operator e~'^"*j4e'^"* are required, and it is this crucial ingredient which allows us 
to obtain, similar to the calculation in [14], the exact result 



AT 




(7) 




(8) 



N/2 



{A){t) = (A)(0) cos{2ht) Y[ cos^[2AA,e(j)t]. 



(9) 
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Fi gur6 1. Time evolution of the expectation value (9) of an observable A for magnetic field h — Q and various system sizes A^. 

Left: For ct — 2, an apparent decay is observed, superimposed by oscillations. The N dependence of the time evolution 
is so weak that the curves for the various system sizes cannot be discerned in the plot. The behaviour for other a > 1 
is qualitatively similar. Right: For a — 1/2, the expectation value again appears to be decaying, but on a time scale 
that depends strongly on the system size N (note the logarithmic scale). Similar behaviour is observed for other values 
of a. between zero and one. 



In comparison with the original Emch-Radin model, the important difTerence here is the explicit dependence 
of the argument of the cosine through the normalization Na- Regarding the approach to equilibrium, the Larmor 
precession cos{2ht) is not relevant and we set h — in the following. The behaviour of (9) is plotted for exponents 
a = 2 and q = 1/2 and various system sizes A'^ in Fig. 1, and in all cases the expectation value of A appears to 
be decaying in time to the microcanonical ensemble average {yl)mic ~ 0. However, this decay is only apparent, 
as we can read off from (9) that {A}{t) is an almost periodic function in time for all finite A*", and recurrences 
(or Loschmidt echos) will therefore occur on much longer time scales than shown. To possibly observe true 
equilibration, we have to invoke, as often in statistical physics, the idealizing concept of the thermodynamic limit. 



4. Thermodynamic hmit 



In the thermodynamic limit, recurrence times may diverge and an relaxation to equilibrium may take place. An 
analysis of the infinite system dynamics is most rigorously done in a *-algebraic language [15, 17], but it essentially 
boils down to discussing the large- A' limit of the product in (9), which is the content of the present section. To 
keep the presentation simple, the following derivations are done for algebraically decaying interaction strengths, 
e(j) = j^" with some positive exponent a, but generalizations are possible. 

4.1. Upper bound on {A){t) 

The aim of this section is to construct, in the thermodynamic limit N ^ oo, an upper bound on the product 



N/2 , . r X 



(10) 
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occurring in the expression for the expectation value {A){t) in (9). Since j cos(x)| < 1, we have 



N/2 

Pn<11 

Then, for given t and A'^ there exists a finite integer 

Mjv(t) = 



1/q 



defined such that 



r ^ 2 



Vj > Mjv(t). 



Hence we can split the product in (11) in the following way, 



A/„{t)-l 

p.v< n 



f2Nct 



N/2 

n 



< 1 



(11) 



(12) 



(13) 



(14) 



For finite and some large t, it will happen that the second product in (14) consists of no factors at all, resulting 
in a trivial (and fairly useless) upper bound Pn < 1. In fact, this was to be expected, since the unitary, quasi- 
periodic dynamics of finite systems gives rise to recurrences. But we are interested in the large- A*' limit where, 
for a given t and large enough A'', the second product in (14) will consist of a large number of factors. By virtue 
of the definition of Afjv(t), the arguments of the cosines in these factors are all between —n/2 and n/2, and for 
this range of arguments the cosine can be bounded by 



I cos(a;)| < 1 



4x^ 



(15) 



From elementary properties of the exponential function it then follows that 



1 - 



n/2 

lim Pn < lim TT 

/V— ^oo iV— >oo 

j = A/„(t) L 

N/2 

< lim TT exp 

iV-KX) 

j=M„(t) 



lim exp 

iV— >oo 



Wat 



2 N/2 



T - 



(16) 



Asymptotically for large A', the sum in (16) is given by^ 



N/2 N/2 

Z. ^ Ala /_ ^ AIol 



for0<a<l/2, 
C(2a) forl/2<a<l. 



(17) 



^ In equations (17) and (18) we exclude for simplicity the case a = 1/2 which, as in (6), would lead to corrections 
logarithmic in N . 
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where the asymptotic formula (6) has been used in the last step. Inserting this result into (16), a Gaussian in t 
bound on Pn is obtained which translates, via (9), into a Gaussian upper bound on the expectation values of the 
observables A, 



f 

exp 




exp 





for < a < 1/2, 
for 1/2 < a < 1. 



(18) 



In both cases, and therefore for all < a < 1, the Ai'-dependent exponent goes to zero in the large-A'^ limit. 
Hence, for long-range interactions, the upper bound is constant in time, limjv^oo j(^)(t)| < |(^)(0)|, and no sign 
of relaxation to equilibrium can be inferred from the bound. In a way, this is not too surprising, as this finding is 
in agreement with the large- A'' trend of the finite-system results plotted in Figure 1. A more detailed discussion 
of the bound (18) will follow in Section 4.3. 

4.2. Lower bound on 

The finite-system results in Figure 1 indicate that, for long-range interactions with exponents < a < 1 and in 
the thermodynamic limit N ^ oo, the Emch-Radin model apparently fails to approach equilibrium. To prove 
that this is indeed the case, and that the failure to construct a nontrivial upper bound in Section 4.1 was not 
due to the crudeness of the bound, we will derive a lower bound on the product Pjv for the case of long-range 
interactions. 

Starting from definition (10), we can write 

N/2 



cos^ 



N/2 

= exp<; ^In 



For any given t there exists an A'o(f) such that 



< 1 



\fN > iVo(t). 



(19) 



(20) 



Then we can use the inequality cos'^(a;) > 1 — a;'^ to obtain the bound 



fN/2 

Pn > exp < In 
lj=i 



1 



2jVct 



valid for A'' > No{t). In the thermodynamic limit we can replace the sum by an integral, 
lim Pjv > lim exp< 



lim exp 

iV— >oo I 



<^ 2Nxa 



(21) 



1 F ( ^ M Kf^^^i"' 
^"'^^'"2^'^'^"2^'liV^ 



2A/"c.f 



=1/2 
J x=1/N 



(22) 
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where 2F1 denotes the Gauss hypergeometric function. The primitive in the second line of (22) was obtained 
with Mathematica. Using the large-A'' asymptotic equaUty 

2Mo.t (l-Q)2^-"t 

iV- ~ N ^ 

as obtained from (6) for < q < 1, the argument of the exponential function in (22) can be rewritten as 



2a 2a' 4"-iiV2-2" y \ 2a 2a m 

^2,2\ / n „\2,2 



2a H- iVa -f — In 1 - — - In 1 - ; ^ ' „ 



(24) 



= 0{1/N) 

From the Gauss series expansion (15.2.1 in [20]) we obtain the large- A?^ expansions 



2«2F. (1, -i^; 1 - ±; = 2a + O (iV^-) , (25) 

Na .F, (1, 1 - ^; ^^^^^) -Na + O (iV^) , (26) 

finding that the leading order of the overall exponent (24) is O {N^^) or O (A'^'^"^^), depending on which one is 
larger. This results in a lower bound 

lim Pn> 1, (27) 
valid for exponents < a < 1. Together with the upper bound (18), this proves that 

hm {A){t) = {^)(0). (28) 

So indeed, as suggested by extrapolating the finite-system plots in Figure 1, the infinite system is trapped for all 
times in its initial state, showing no sign of relaxation to equilibrium. Having in mind physical systems of finite 
size, it seems more useful to state the result in the following way. 

Proposition 4.1. 

Consider the Ising-type Hamiltonian Hn defined in (2) with power law decaying interactions e{j) = j~" with 
< a < 1. Consider further an observable A of the type (7) and an initial state p(0) bemg diagonal in the 
af tensor product eigenbasis of the underlying Hilbert space H. Then for the expectation value {A){t) of A with 
respect to pit), the following holds true: For any fixed time r and some small 5 > Q, there is a finite no(r) such 
that 

\{A){t)-{A){Q)\<5 \ft<r, N >no{T). (29) 



Interpreting this result in terms of an experiment, we can think of an experimental resolution S for the measurement 
of A, and some duration r of the experiment. Then the above proposition states that, within the experimental 
resolution and for a large enough system, no deviation of {A){t) from its initial value can be observed for times 
t < T. 
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Figure 2. Time evolution of the expectation value (9), plotted as a function of (the logarithm of) for various system sizes N. 

The graphs for different are hardly discernable in the plot, as they collapse onto each other. The three plots are for 
different values of the exponent a (from left to right: a — 0, 1/4, 3/4). The dashed lines show the upper bounds as 
given by Equation (18). 



4.3. Relaxation to equilibrium in rescaled time 

Although for the long-range Emch-Radin model in the thermodynamic limit relaxation to equilibrium does not 
take place for any finite time t, the plots in Figure 1 suggest that equilibration can be observed when considering 
a suitably rescaled time variable. The case of Curie- Weiss-type interactions, i.e. an exponent a — 0, can serve as 
a simple, exactly solvable example: Defining the rescaled time variable tq = t/^/N , the product Pn in (fO) can 
be written as 



exp < A'^ In 



2ro 



= exp{-4ro' + 0(f/iV)} 



(30) 



yielding a Gaussian relaxation to equilibrium in the large- A limit. For general exponents < a < 1, such a 
simple calculation is not feasible. However, comparing the upper bound (18) with the corrections to the lower 
bound in (25) and (26), we can read off that an a-dependent rescaled time 



tA-i/^ for0<a<l/2, 
tN°'-^ for 1/2 < a < 1, 



(31) 



leads to A-independent Gaussian bounds. This implies that the expectation value I^A)(tc,) displays a Gaussian 
decay to equilibrium when plotted as a function of r^, and the microcanonical expectation value is correctly 
reproduced in the limit of — >■ oo, 

lim lim (A)(r„) = 0. (32) 

We checked this prediction by numerically computing the product Pn in (10) for various values of A, and plotting 
the result as a function of r^. As is shown in Figure 2, the graphs of Pn for different A nicely collapse onto a 
single curve when plotted versus rescaled time Tq. For comparison, the bounds (18) are plotted as dashed lines 
in the same figure. This confirms that the time scale predicted by (31) is indeed the correct one. Rephrasing this 
result, we have found that the lifetime To of quasi-stationary states in the long-range Emch-Radin model diverges 
with the system size like To ~ A'^ with an exponent q — min{l/2, 1 — a}. 
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5. Discussion and outlook 

We have analytically studied the time evolution of long-range quantum spin models with Ising-type interactions 
(2) on one-dimensional lattices and for a certain class of initial states. Generalizations to higher-dimensional 
lattices are straightforward, but have not been discussed in this article. For finite system sizes A^, almost periodic 
behaviour is observed with recurrence times that increase exponentially with TV. In a rigorous sense, equilibration 
can therefore occur only in the thermodynamic limit of infinite lattice size. Whether this indeed happens depends 
on the asymptotic behaviour of the interaction strength e{j) for large distances j on the lattice. To analyze the 
long-time asymptotic behaviour of the expectation value {A){t) in the thermodynamic limit, we constructed upper 
and lower bounds. For short-range interactions, i.e. interactions decaying faster than Cj~" with some constant 
C and exponent q > 1, it had already been proven in [15] that relaxation to equilibrium takes place as expected 
and the expectation value {A){t) approaches the microcanonical ensemble average (yl)mic — for large t. 
For exponents a between zero and one, the interaction strength is not absolutely summable and an A''- 
dependent normalization Ma is needed in (2) to render the energy per spin finite. In this case, the trivial (i.e. 
constant-in-time) upper bound {A){t) < {A){0) is obtained, which does not permit any conclusions on whether 
or not (A) approaches its equilibrium value (^)mic ~ 0. In order to get a more detailed picture of what happens 
in this case, we constructed a lower bound, finding that {A){t) > (^)(0) for exponents < a < 1. Rephrasing 
this result, our Proposition 4.1 asserts that, at a given instant of time and for large enough A'', {A){t) will be 
practically unchanged from its initial value (^)(0). This means that, for large enough systems, equilibration 
occurs on a time scale beyond the experimentally accessible one and will not be observed in practice. Despite 
the superficial similarities, such a behaviour is notably different from the failure of thermalization reported for 
near-integrable systems [21-23]. Our results extend the concept of quasistationary states, previously observed 
and extensively studied in classical systems, into the realm of quantum mechanics. Moreover, we provide the 
first analytic proof of quasistationary behaviour from the microscopic time-evolution equations, as well as an 
expression for the A'^-dependence of the corresponding lifetime. To ~ Af' with q — min{l/2, 1 — a}. 
It is tempting to speculate that quasistationary behaviour, i.e. equilibration times that diverge for large system 
sizes A'', might show up under much more general conditions, or even generically for long-range interacting quantum 
systems. An extension of our results to general anisotropic Heisenberg models might be a promising future project 
and, although an exact solution for the time evolution of (A) seems to be out of reach, the derivation of suitable 
bounds on (A) (t) might be feasible. Finally, the simultaneous presence of short- and long-range interactions might 
lead to equilibration taking place on two different time scales, possibly leading to more complex behaviour. The 
peculiarities of diverging relaxation times are expected to be reflected also in many applied aspects of quantum 
spin systems, including the study of quantum quenches [24, 25], decoherence, and many others. 
As a concluding remark, a comment on the role of the normalization factor A/'a in the Hamiltonian (2) is in order. 
Since the Hamiltonian enters the time-evolution equation, it is not surprising to see that the system's dynamics 
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inherits the A^-dependence of Ma- But this is not the only source, as the term discussed in (17) introduces a 
further A*'-dependent contribution, and it is the interplay of the two terms which determines the nontrivial TV- 
scaling of the relaxation time. So how should the prefactor Ma in the Hamiltonian then be chosen? The answer to 
this question depends on the physical problem under investigation, and, depending on the situation, the following 
three choices may have their virtues: 

1. The normalization factor Ma as defined in (5) renders the energy an extensive quantity. As a result, 
energy and entropy have the same A^-scaling in the large- limit. This is a necessary requirement for 
having competition between energetic and entropic effects, which in turn is crucial for the occurrence of 
thermodynamic phase transitions. 

2. An A'^-dependent prefactor which leads to A'^-independent long-time asymptotic dynamics. This would 
eliminate quasi-stationary behaviour, but it is a very impractical choice, as it would require knowledge of 
the long-time asymptotic behaviour. 

3. A prefactor independent of A'', as this is what is realized in most physical systems. Such a choice is of course 
the most useful one when comparing experimental and theoretical results for systems of a given finite size, 
but it has the disadvantage of being unsuitable for the calculation of thermodynamic limits. 
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